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Eric Braaten 



D i The short-distance coefficients for the color-octet term in the fragmenta- 



tion function for a gluon to split into heavy quarkonium states is calculated 
^ \ to order a^. The gauge-invariant definition of the fragmentation function by 

Collins and Soper is employed. Ultraviolet divergences are removed using 
the MS renormalization procedure. The longitudinal term in the fragmen- 
tation function agrees with a previous calculation by Beneke and Rothstein. 
The next-to-leading order correction to the transverse term disagrees with a 
previous calculation. 
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Typeset using REVTgX 



I. INTRODUCTION 



The cross sections for heavy quarkonium states probe the production of heavy-quark- 
antiquark pairs with small relative momenta. Many of the theoretical uncertainties in 
quarkonium production decrease at large transverse momentum. Factorization theorems 
for inclusive single-hadron production ||T[] guarantee that the dominant mechanism for pro- 
ducing heavy quarkonia with high px is fragmentation 0], the production of a parton which 
subsequently decays into the quarkonium state and other partons. This process is described 
by a fragmentation function D{z,fi), where z is the longitudinal momentum fraction of the 
quarkonium state and /i is a factorization scale. 

The NRQCD factorization formalism |^ can be used to factor the fragmentation func- 
tions D{z,n) for quarkonium into NRQCD matrix elements, which can be regarded as phe- 
nomenological parameters, and short-distance factors, which depend on z and are calculable 
in perturbation theory. Most of the phenomenologically relevant short-distance factors begin 
at order and have been calculated to leading order. However there is one short-distance 
factor that begins at order a^. It is the color-octet ^5*1 term in the gluon fragmentation 
function, whose NRQCD matrix element is denoted (Osi'^Si)). This term is of particu- 
lar phenomenological importance. Braaten and Yuan showed that it must be included in 
the gluon fragmentation function for triplet P-wave states in order to avoid an infrared di- 
vergence in the short-distance coefficient of the color-singlet matrix element {Oi{^Pj)) 0. 
Braaten and Fleming argued that the {Os{^Si)) term is also phenomenologically necessary 
in the gluon fragmentation function for spin-triplet S-wave states in order to explain the 
production rate of direct J/ip and tp' at large pt at the Tevatron This led to the re- 
markable prediction by Cho and Wise that J/ip and ip' at large pt should be transversely 
polarized 0. 

In the earliest calculations of fragmentation functions for heavy quarkonium the 
short- distance factors were deduced by comparing the cross sections for quarkonium pro- 
duction with the form predicted by the factorization theorems for inclusive single-hadron 
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production. However the fragmentation functions can also be defined formally as matrix el- 
ements of bilocal operators in a light-cone gauge |^ or, more generally, as matrix elements of 
non-local gauge- invariant operators The gauge- invariant definition of Collins and Soper 
was first applied to calculations of the fragmentation functions for heavy quarkonium by Ma 
TI| . The definition is particularly convenient for carrying out calculations beyond leading 



order in as- It was used by Ma to calculate the short- distance factor of the color-octet ^5*1 



term in the gluon fragmentation function to next-to-leading order in |Tl 



In this paper, we calculate the short- distance coefficients for the color-octet "^Si term in 
the fragmentation function for a gluon to split into heavy quarkonium states to order a^. 
We use the gauge-invariant definition of the fragmentation function given by Collins and 
Soper [||], and we remove ultraviolet divergences using the MS renormalization procedure. 
Our result for the longitudinal term in the fragmentation function agrees with a previous 



calculation by Beneke and Rothstein [12|. Our result for the next-to- leading order correction 



in the transverse term disagrees with a previous calculation by Ma |11 



II. GAUGE-INVARIANT DEFINITION 

The fragmentation function Dg^H^z, ^) gives the probability that a gluon produced in 
a hard-scattering process involving momentum transfer of order /i decays into a hadron H 
carrying a fraction z of the gluon's longitudinal momentum. This function can be defined 
in terms of the matrix element of a bilocal operator involving two gluon field strengths in a 
light-cone gauge 0. In Ref. Collins and Soper introduced a gauge- invariant definition 
of the gluon fragmentation function that involves the matrix element of a nonlocal operator 
consisting of two gluon field strengths and eikonal operators. One advantage of this definition 
is that it avoids subtleties associated with products of singular distributions. The gauge- 
invariant definition is also advantageous for explicit perturbative calculations, because it 
allows the calculation of radiative corrections to be simplified by using Feynman gauge. 

The gauge-invariant definition of Collins and Soper is 
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FIG. 1. Leading order Feynman diagram for g QQ. 



N-2 



dx e 



IQtx{N - l)k+ 



(1) 



The operator £{x ) in ([I|) is an eikonal operator that involves a path-ordered exponential 
of gluon field operators along a light-like path: 



£{x- 



Iba 



P exp 



+ig / dz^A-^{0+,z-,0^) 



ba 



(2) 



where A^{x) is the matrix-valued gluon field in the adjoint representation: [A'^(x)]ac = 
The operator 'Ph{p+,pj^) in (13) is a projection onto states that in the asymptotic 
future contain a hadron H with momentum p = {p^,p~ = {mfj + p'j_)/p'^ ,p±), where ttih is 
the mass of the hadron. In the definition (|l|), the hard-scattering scale fi can be identified 
with the renormalization scale of the nonlocal operator. In perturbative calculations, it 
is convenient to use dimensional regularization to regularize ultraviolet divergences. The 
prefactor in the definition (|l]) has therefore been expressed as a function of the number of 
spatial dimensions = 3 — 2e. 

For any state H that can be defined in perturbation theory, the definition (|I|) can be used 
to calculate the fragmentation function Dg^niz, fi) as a power series in a^. A convenient 
set of Feynman rules for this perturbative expansion is given in Ref. [Q]. If the state consists 
of a QQ pair with invariant mass p^, the lowest order diagram is shown in Fig. |l]. The 



circles connected by the double pair of lines represent the nonlocal operator consisting of 
the gluon field strengths and the eikonal operators. The momentum k = (fc"*", k~ , k±) flows 
into the circle on the left and out the circle on the right. The cutting line represents the 
projection onto states that in the asymptotic future include a QQ pair with total momentum 
p = {zk~^,p'^/ {zk~^), 0_l). If the heavy quark has relative momentum q in the QQ rest frame, 
the invariant mass is p"^ = 4:{mQ + q^). 

The fragmentation of a gluon into heavy quarkonium state H involves many momentum 
scales, ranging from the hard-scattering scale /i, which we will assume to be larger than tuq, 
to momenta much smaller than mq where nonperturbative effects are large. The NRQCD 
factorization formalism allows the systematic separation of momentum scales of order mq 
and larger from scales of order mqv or smaller, where v is the typical relative velocity of the 
heavy quark in the hadron. The factorization formula has the form 

= Y.dUz^t^){^nJ ■ (3) 

mn 

The short- distance coefficients dmn{z,^) are independent of the quarkonium state H and 
can be calculated as a perturbation series in as{mq). All long-distance effects are factored 
into the NRQCD matrix elements (O^^), which can be expressed as matrix elements in an 
effective field theory. They have the general form 

{OD = (0|x^/C> Vh ^^IC^m , (4) 

where JCm and /C„ are constructed out of color matrices, spin matrices, and covariant deriva- 
tives and the operator Vh projects onto states that in the asymptotic future contain a 
quarkonium state H. The NRQCD matrix elements are nonperturbative but they are uni- 
versal, with the same matrix elements describing inclusive production in other high energy 
processes. 

The threshold expansion method |T3|JT^ is a general prescription for determining the 
short- distance coefficients dnm{z.,^) in the NRQCD factorization formula (^. The diagrams 
for the fragmentation function of a QQ pair are computed in perturbation theory, except 
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that the QQ pair is allowed to be in a different state on the two sides of the final-state cut. 
To the left of the cutting line, the Q and Q have relative momentum q in the QQ rest frame 
and color and spin states specified by Pauli spinors rj and ^. To the right of the cutting line, 
the Q and Q have relative momentum q' and color and spin states specified by rj' and 
After expanding around the threshold q = q' = 0, the resulting expression for the diagrams 
has the form 

(5) 

mn 

where the color and spin matrices and are polynomials in the relative momenta q and 
q'. The short-distance coefficients dnm{z, ^i) in the factorization formula @) can then be read 
off from this expression. For example, if the sum of the diagrams includes the color-octet 
^Si terms 



Nrric 



dT{z){6'^ - z'z^) + dL{z)z'z^ i^a^T^i' ^^a'T'^r], (6) 



then the fragmentation function includes the terms 
N 



Dg^H{z) 



Arar 



drizW^ - z'z^) + dL{z)z'z^ (0|xV^T>Ph V'^fxT'^xlO). (7) 



If we sum over the spin states of H, the matrix element in (|^) is proportional to 5^^ and (|^ 
reduces to 



D,^h{z) = [{N - l)dT{z) + dL{z)] {O^i'S,)), (8) 

where {Oq{^Si)) is the color-octet ^5*1 matrix element defined in Ref. 0: 

{O^CS,)) = ^(0|xV^T>P^V'V^T'^x|0). (9) 

The factor of l/(4mQ) accounts for the relativistic normalization of the projection operator 
used in Refs. [111111. 
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III. LEADING ORDER 



The only Feynman diagram of order for the fragmentation process g QQ is shown 
in Fig. |l]. Using the Feynman rules of Ref. [|l|, the expression for the diagram can be 
easily written down in terms of spinors u and v that describe the color, spin, and relative 
momentum states of the Q and Q. We can allow for the Q and Q on the right side of the 
final-state cut to have different color and spin states and different relative momentum q' by 
replacing their spinors by u' and v' . The resulting expression is 



Setting q = q' = 0, we can replace the Dirac spinors with Pauli spinors by substituting 

u-iJ^'v = 2mQ L^i eV'T^, (11) 
where L^i is a boost matrix that satisfies the identities 

- = (12) 

n°L„, = (13) 
VP 

The expression (|T0|) then reduces to 

Si^-^IWQ ~ ~ z^z^V^^'Tr eV^T^. (14) 

Comparing with (P), we can read off the order-a^ terms in the short-distance functions dxiz) 
and (Il^z) defined in (0): 

d'!:°Hz) = 0. (16) 

The dependence on the number of spatial dimensions agrees with that in Ref. il^ . 
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(b) 






(d) J (e) J 

FIG. 2. The Feynman diagrams of order for g QQ with QQ final states. There are 
additional contributions from the complex-conjugate diagrams. 

IV. VIRTUAL CORRECTIONS 



The Feynman diagrams for the fragmentation function for g QQ at order a1 consist 
of virtual corrections, for which the final state is QQ, and real-gluon corrections, for which 
the final state is QQg. The diagrams with virtual-gluon corrections to the left of the cutting 
line are shown in Fig. H. The black blob in Fig. |^(a) includes the vertex corrections and 
propagator corrections shown in Fig. ^. 

We calculate the diagrams using Feynman gauge. In this case, the diagram Fig. ^(b) 
vanishes, because the gluon attaching to the eikonal line gives a factor of n^. When con- 
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(a) 





(c) 



FIG. 3. One loop correction diagrams for g* — > QQ. 

tracted with the factor k ■ ng'^^ — p^n^ from the circle on the left side of the cut, it gives 
a factor {k ~ p) ■ n = 0. In the limit q = q' = 0, all the other diagrams can be reduced 
to the leading order diagram in Fig. |l| times a multiplicative factor. For the diagrams in 
Fig. |](a), this follows from the Dirac equation for the heavy-quark spinors. For the remain- 
ing diagrams, we must also use the fact that the contraction of n'^ with the right side of the 
diagram vanishes. 

The virtual corrections contribute only to the transverse short- distance function d^iz) 
defined in (0). We will express the various contributions in the form of the leading-order 
result ([TB|) times a multiphcative factor. The sum of the diagrams in Fig. 0(a), together 
with their complex conjugates, is 



S^''\z) = d)r""{z) X 2 Re 



(LO), 



A + n+(Zo-ll 



(17) 



where A is the vertex correction factor from the diagrams in Fig. 0(b), 11 is the propagator 
correction factor for a gluon with invariant mass 4mQ from the diagrams in Fig. |^(d), and 



Zq — 1 comes from the wavefunction renormalization factors Z]^"^ for the heavy quark from 
the diagrams in Fig. ^(c). These correction factors are given by 



A 



2\ ^ 



71 \m%, 



12 evv 12 €;r 6 2 

_l£(l±i)_2I(i±i)_l_21u2 

3 e/if^ 3 



3 euv 
15 - r(l + e 



+ 



93 - lOn 



(18) 
(19) 
(20) 



12 euv 36 
where nj is the number of hght quark flavors. The subscripts on the poles in e indicate 
whether the divergences are of ultraviolet or infrared origin. The sum of the virtual correc- 
tions given in (0) is 

■2^' no ^_ 3 123 -lOn 



dr)(z) = 4'°^(.)x^ p^) r(l + 6) 



7r V mr 



12 -Uf 



3e 



uv 



2e 



+ 



+ 21n2 



IR 



(21) 



Ma's result for these diagrams |]TT| differs by an additive constant —5.487 inside the square 
brackets and by setting A^(A^ - 1) 3(iV - 1) in d^^^\z). 

The contribution from the diagrams in Fig. |](c) with its complex conjugate is 



42-) (^) = 4^^^(z) X 247ra,/i^^ Re 



7{LO) 



i {2k ■ uIabc - Iab) 



(22) 



where the scalar integrals Iab- - are given in appendix A. The contribution from the sum of 
the diagrams in Figs. @(d) and 0(e) with their complex conjugates is 



(2d,e). 



d^T^\z) X 967iasfiVQ Re 



i {k ■ nl 



ABCD — J^ABD) 



(23) 



The sum of the virtual corrections given in (12^) and (EBI) is 



d. 



(2c,d,e) 



.)=4'"^Wx^f^Vr(l + e) 



3(1 -e) 



uv 



+3 + 61n2 + 61n^2 

2 



(24) 



Ma's result for these diagrams |TT[ differs by an additive constant —3 + ^vr — 6 In 2 inside 
the square brackets and by setting A^(A^ — 1) 3{N — 1) in dj:''^\z). 
The total virtual correction is the sum of f^Tl) and (123): 
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J (e) 

FIG. 4. The Feynman diagrams of order for g QQ with QQg final states. There are a 
total of 25 diagrams, but only the left halves of the diagrams are shown. 



Vr \ 1717 



3(l-6)r(l + e) r(l + e) 



^UV^IR 



+ 



177 - lOn 



8 In 2 + 6 In^ 2 



(25) 



where (3q = (33 — 2nf)/Q. 



V. REAL GLUON CORRECTIONS 



The Feynman diagrams for the real-gluon corrections to the fragmentation function for 
g QQ pair are shown in Fig. ^. We draw the 5 left-half diagrams only, but they must 
be multiplied by their complex conjugates to give a total of 25 diagrams. Only diagrams 
^(a), ^(b) and ||(c) contribute to the color-octet ^Si term, which reduces the total number 
of diagrams to 9. We calculate the diagrams using Feynman gauge. After a considerable 
amount of algebra, they reduce to 
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^'^^"8iV(iV-l)m| 7rr(l-e) IvmlJ 1^' ^(1-^)/ 4-.)// ' ^^^^ 

where t = (1 — + 2; — 1), x = 2g • p/p^, q is the final-state gluon momentum, and p is 

the QQ momentum. Integration over x can be done using the following identities: 

00 1 _ p 

dx^ = z{l - zy-^' - — - B{1 + e, 1 - e) , (28) 

(i-z)/z euv 

dx^— = z(l - z)'^-^' B(l + e, 1 - e) . (29) 

il-z)/z X^ 

The subscript on e in ( pSf ) indicates that the pole has an ultraviolet origin. In (|26|), there 
is also an infrared divergence associated with the limit z 1. It can be made explicit by 
using the expansion: 

The final results for the real-gluon corrections are 

3(1-6 



dp^\z) = , X ^ f^yr(l + e) 

3(1 — e) / z 1 — z . . 

' ^ ^ ' ' +z(l~z) 



6{l~z) 



6 /ln(l 



z 



z \ 1 — z 



+ Q{2~z + z^)\n{l-z 



(31) 



,(real)/ x TTa^ 3as 1 - ^ /„„x 

d} '(z) = 5- X . (32) 



Note that the double-pole term proportional to l/ieuv^m) in ( pl| ) exactly cancels its counter 
part in the virtual correction (pS]). The sum of ( pSj ) and ( |3TD is free of infrared divergences. 
Ma's result for the real-gluon corrections is completely different from the sum of ( |3TD 
and (P^). In particular, he found that the [ln(l — z)/{l — z)]+ terms cancelled. In the way 
we organized the calculation, there is no possibility of such a cancellation. 
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VI. RENORMALIZATION 



The sum of the order-a^ corrections to d^iz) in (^) and (|3T|) still contains ultravio- 
let divergences in the form of single poles in e. These divergences are cancelled by the 
renormalization of the coupling constant as in the leading-order expression ([T6|) and by the 
renormalization of the nonlocal operator in (^. The renormalization of in the MS scheme 
can be carried out by making the following substitution in (|TBp: 

1 



a., a.. 



(33) 



where f3o = (33 — 2n/)/6. The operator renormalization in the MS scheme can be carried 
out by making the following substitution 

1 f^dy 



ZTT 



euv Jz y 



P,,{y) d^^''\z/y) 



(34) 



where Pggiy) is the gluon splitting function: 

l-z 



/9o 



l-z) 



+ + z{l-z) + '^ 5(1 - z) 

Z 



(35) 



The sum of the two contributions of order a1 from renormalization is 



dp''\z) 



2e 



m{N - l)ml 



(-3) 



(36) 



This cancels the single ultraviolet poles in the sum of (p5D and (|3lD. Our final result for the 
transverse fragmentation function is obtained by adding the order-a^ corrections to dxiz) 
from dl), (|1]), and and taking e ^ 0: 

1^ 



driz^n) = ^ < 5(1 - z) + 



48m? 



A(/i)5(l -z)+{\n 



2mQ 2^ 



+Q{2 - z + z^)\yi{1 - z) - 



6 /ln(l 



z \ 1 — z 



(37) 



where the coefficient A{fi) is 
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A{fi) = /3o In 



jj 13 
2mQ 6 ^ 



+ 81n2 + 61n^2 



(38) 



The result ( ]37| ) disagrees with the final result obtained by Ma |Tl|]. Our final result for the 



longitudinal fragmentation function is obtained by setting e ^ in (|32D : 

alifi) 1 - z 



(39) 



Q 



This agrees with the result of Beneke and Rothstein |T2|. The fragmentation probabilities 
obtained by integrating dxiz) and dL{z) diverge, because these functions behave like 1/^; as 
z ^ 0. The higher moments of the fragmentation functions however are well-defined. 



VII. DISCUSSION 

The color-octet ^5*1 term in the fragmentation function is important for calculating the 
production at large pt of spin-singlet S-wave states, like the J/V', and spin-triplet P-wave 
states, like the Xcj- We can deduce the fragmentation functions for each of their spin states 
by using the approximate spin symmetry of NRQCD to simplify the expression (^. The 



functions dxiz) in (|37[) and dL{z) (^) give the fragmentation functions for the transverse 
and longitudinal spin states of the J/ip, respectively: 

/^5-^/^(±i)(^) = dT{z){Oi'\'S,)) , (40) 
D,^jl^^,){z) = d,{z){Oi/\'S,)). (41) 

The sum over spin states is 

D,^j,.^{z) = [2dT{z) + dL{z)] {Oi/'^i'S^)) . (42) 



The fragmentation functions for each of the spin states of the Xcj are [16 



D,^^^,{z) = [2dT{z) + d^z)] ion's,)) , (43) 

D,^^^,^o){z) = 3dT{z){OreS,)), (44) 

= I [dAz) + dLiz)] (ores,)) , (45) 
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^™(o)(^) = [driz) + 2dL{z)] {OrCS^)) , (46) 
= I [driz) + dr.{z)] {OreS^)) , (47) 
^™(±2)(^) = 3dT(^)(08^=°(=^^i)) . (48) 

The sums over spin states are 

D,^^^^{z) = (2J + 1) [2dT{z) + dL{z)] {OreS^)) . (49) 

In order to give accurate predictions for the production of quarkonium at large pr, it is 
important to know the next-to-leading order correction to the color-octet ^Si term in the 
gluon fragmentation function. Unfortunately our calculation of the short-distance coefficient 
disagrees with the previous calculation by Ma. An independent calculation of this important 
function is therefore essential. 

This work was supported in part by the U.S. Department of Energy Division of High 
Energy Physics under grant DE-FG02-91-ER40690, by the Alexander von Humboldt Foun- 
dation, and by the Korea Institute for Advanced Study. J.L. would hke to thank the OSU 
theory group for its hospitality during his stay in Columbus. 



APPENDIX A: INTEGRAL TABLE 

In this appendix, we present the explicit values of the integrals encountered in evaluating 
the virtual-gluon corrections. These integrals have the form 

f d^^H 1 

= i (27r)^+i AB... ' ^^^^ 
where the denominator AB ■ ■ ■ can be a product of 1, 2, 3, or 4 of the following factors: 

A = f + ie, (A2) 

B ^{l-pf + ie^l'^ -21-P + 4m| + ie, (A3) 

C = (p-l) -n + ie, (A4) 

D ^ {l-p/2f -ml + ie^f -l-p + ie. (A5) 
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The momentum p is that of a QQ pair with zero relative momentum (jp — Auiq) and n 
is Ught-Uke (n^ = 0). The integrals I a and Ib vanish in dimensional regularization. By 
symmetry under p — > / — p, we have Iad = Ibd- Some of the integrals can be reduced to 
ones with fewer denominators by using the identity A-\- B — 2D — Attiq: 



AuiqIabd = 2{Iad - I as), 



AmnlABCD — IacD + I BCD — '^^ABC- 



The independent integrals that need to be evaluated are therefore 
i /vre'^V r(l + e)r2(l -e) 



Iab 
Iad 
Iabc 
Iacd 
Ibcd 



(47r)2 \^ ml 
i f Att 



ef7yr(2 - 2e) 

r(i + e) 



{ATry WqJ euvil-2e) 



—I 



Tie 



r(i + e)r2(i 



(47r)2p -nyml J euveiRT{l - 2e) 
+i /47rV r(l + e) 



(47r)2p ■ n \ml 



TT 



21n2 + e ( y -61n^2j +0(e^) 



47rY r(l + e) 



(A6) 
(AT) 



(A8) 
(A9) 
(AlO) 
(All) 
(A12) 



{AttYp ■ n \mQ ) euv^m 

The subscripts on the poles in e indicate whether the divergences are of ultraviolet or infrared 
origin. 
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